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Abstract 

A search for the lepton flavour violating processes ep — > fiX and ep — > tX is performed 
with the HI experiment at HERA. Final states with a muon or tau and a hadronic jet are 
searched for in a data sample corresponding to an integrated luminosity of 66.5 pb" 1 for 
e + p collisions and 13.7 pb _1 for e~p collisions at a centre-of-mass energy of 319 GeV. No 
evidence for lepton flavour violation is found. Limits are derived on the mass and the cou- 
plings of leptoquarks inducing lepton flavour violation in an extension of the Buchmuller- 
Riickl-Wyler effective model. Leptoquarks produced in ep collisions with a coupling strength 
of A = 0.3 and decaying with the same coupling strength to a muon-quark pair or a tau- 
quark pair are excluded at 95% confidence level up to masses of 459 GeV and 379 GeV, 
respectively. 
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1 Introduction 



In the Standard Model (SM) the particle interactions conserve lepton flavour, although there 
is no underlying symmetry supporting this feature. However, experimental evidence for lepton 
flavour violation (LFV) in solar and atmospheric neutrino oscillations has been reported [1,2]. 
The experimental upper bounds on neutrino masses imply very small LFV effects in the charged 
lepton sector. The observation of such effects would clearly indicate new phenomena beyond 
the SM. 

In ep collisions at HERA, LFV processes ep — > fiX or ep — > rX lead to final states with a 
muon or a tau and a hadronic system X. The LFV process can proceed via the exchange of a 
leptoquark (LQ), a boson with both lepton and baryon quantum number which appears naturally 
as a colour triplet scalar or vector boson in many extensions of the SM such as grand unified 
theories [3], supersymmetry [4], compositeness [5] and technicolor [6]. 

In this paper a search for LFV phenomena is performed in ep collision data recorded dur- 
ing the years 1998-2000 by the HI experiment, corresponding to an integrated luminosity 
of 66.5 pb" 1 for e + p collisions and 13.7 pb _1 for e~p collisions at a centre-of-mass energy 
y/s = 319 GeV. The present results supercede those derived in previous searches at the HI 
experiment using e + p collisions at s/s = 300 GeV [7]. 



2 Experimental conditions 

A detailed description of the HI detector can be found in [8]. In the following, only the detector 
components relevant for this analysis are briefly discussed. The origin of the HI coordinate 
system is the nominal ep interaction point, with the direction of the proton beam defining the 
positive z-axis (forward direction). Transverse momenta and azimuthal angles are measured in 
the xy plane. The pseudorapidity is related to the polar angle 9 by 77 = — lntan(#/2). 

A tracking system consisting of central and forward detectors is used to measure charged 
particle trajectories and to determine the interaction vertex. The central tracker is composed 
of two concentric cylindrical drift chambers providing full acceptance for particles in the range 
22° < 9 < 160°, complemented by a silicon vertex detector [9] covering the range 30° < 9 < 150°. 
Transverse momenta (P T ) are determined in the central region from the curvature of the particle 
trajectories in a magnetic field of 1.15 T with an effective resolution of 
<j(Pt) I ' Pt — 0.01 • Pr(GeV). The tracking is complemented in the forward region 7° <9< 25° 
by a system of drift chambers perpendicular to the beam axis. 

With a polar coverage of 4° < 9 < 154° and full azimuthal acceptance, the liquid Argon (LAr) 
calorimeter encloses the tracking chambers. It consists of an inner electromagnetic part with 
a fine granularity and an outer hadronic part with a coarser granularity. The energy resolution 
of the LAr calorimeter for electrons and hadrons was determined in test beam measurements to 
be a/E = 12%/v/£(GeV) © 1% and a/E = 50%/^E(GeV) © 2%, respectively [10]. In the 
backward region 153° < 9 < 178°, the LAr calorimeter is complemented by a lead-scintillating 
fibre spaghetti calorimeter. 
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The iron return yoke of the magnet is instrumented with streamer tubes to identify muon 
tracks. Further chambers of the central muon system are positioned around the yoke to provide 
a precise muon track measurement in the polar range 5° < 9 < 175°. Additional drift chambers 
positioned at either side of a toroidal magnet are employed to detect muons in the forward 
direction (3°<#<17°). 

The luminosity is determined from the rate of the Bethe-Heitler process ep—^epj, measured 
using a photon detector located close to the beam pipe at z = —103 m. 

Electrons are identified as compact and isolated electromagnetic clusters in the calorime- 
ter. Within the acceptance of the tracking detectors, an associated track is required. A muon 
candidate is identified by associating an isolated track in the forward muon system or in the 
inner tracking system with a track segment or an energy deposit in the instrumented iron. 
The hadronic final state is reconstructed from the deposits in the LAr calorimeter in combi- 
nation with tracking information. The hadrons are then combined into jets using the inclusive 
£/r-algorithm [11] with a P T -weighted recombination scheme where jets are treated as massless 
and the separation parameter is set to one. 



3 LFV phenomenology and SM background processes 

The LFV processes ep — > fiX and ep — > rX can be attributed to LQs produced at HERA pre- 
dominantly by electron-quark fusion. 

In the framework of the Buchmuller-Ruckl-Wyler (BRW) effective model [12], LQs are 
classified into 14 types with respect to the quantum numbers spin, isospin and chirality. Lepto- 
quarks carry both lepton (L) and baryon (B) quantum numbers. The fermion number F = L+3B 
is assumed to be conserved, taking values of F = 2 for e~q processes and F = for e + q pro- 
cesses. Leptoquark processes proceed via s channel resonant LQ production or u channel virtual 
LQ exchange, as shown in figure 1. For LQ masses m^Q well below the e ± p centre-of-mass 
energy, the s channel production of F = 2 (F = 0) LQs in e~p (e + p) collisions dominates. For 
LQ masses above 319 GeV the s and u channel processes become of equal importance and both 
e~p and e + p collisions have similar sensitivity to virtual effects from F = 2 LQs as well as from 
F = 0LQs. 

The BRW model assumes lepton flavour conservation (LFC) such that the LQs produced in 
ep collisions decay only to eX or v e X final states. These LQs are referred to in the following 
as first generation LQs and have been studied in a recent HI publication [13]. A general ex- 
tension of the BRW model allows for the decay of LQs to final states containing a lepton of a 
different flavour, i.e. p, or r, and a jet, as illustrated in figure 1. Non-zero couplings X egi to an 
electron-quark pair and \ mj (X Tq ) to a muon(tau)-quark pair are assumed. The indices 1 i and j 
represent quark generation indices, such that \ eQi denotes the coupling of an electron to a quark 
of generation i, and X iq is the coupling of the outgoing lepton £ to a quark of generation j. 

'in the following the quark generation indices are attached only when it is relevant. 
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Figure 1 : Left: s-channel resonant LQ production and decay to a lepton-quark pair. Right: 
u-channel exchange of a LQ. The indices i and j represent quark generation indices, such that 
X eqi denotes the coupling of an electron to a quark of generation i, and \t qj is the coupling of 
the outgoing lepton £ to a quark of generation j. For I — pL,T, the LQ introduces LFV. 



The double differential cross section for the s-channel tree level process is [12]: 

fa^ = J_ X 2 eq X 2 £q s 2 [\ scalar LQ 

dxdy 32tts js* - mj Q y + [2(1 - y) 2 vector LQ, K} 

phase space Breit-Wigner LQ propagator term parton density 

where x is the Bj0rken scaling variable, y denotes the inelasticity of the ep scattering process, 
s = sx represents the square of the eq centre-of-mass energy and Flq is the total LQ width. A 
similar expression holds for the u channel exchange [12]. 

An overview of the extended effective model for the LQ coupling to u and d quarks is given 
in table 1. For convenience only one LFV transition is considered: either between the first 
and the second generations or between the first and the third generations. The branching ratio 

LQ — > fx(r)q is given by 

BR = ft x (3 LFV with (3 LFV = - _ and IT, = m LQ X% x { * ^ (2) 

1 rtr) q + 1 e vector LQ 

where Y^ q denotes the partial LQ decay width to a lepton £ = e, fx,r and a quark q and where 
ft = Yf> q l (T £q + T Uiq ) is the fraction of decays into charged leptons. Some LQs, namely Sq, S[, 
Vq and Vi, can decay to a neutrino-quark pair resulting in ft = 0.5. Since neutrino flavours 
cannot be distinguished with the HI experiment, such final states are not covered in this search, 
but they are implicitly included in the search for first generation LQs [13]. 

To determine the signal detection efficiencies, events with LQs are generated using the 
LEGO [14] event generator with the CTEQ5L parametrisation of the parton distribution func- 
tions (PDF) of the proton [15]. The LQ signal expectation is a function of the LQ type, mass, 
coupling constant and /3lfv- The analysis usually requires a large number of simulated signal 
Monte Carlo (MC) samples. To overcome this technical difficulty, the LEGO program is used 
to produce a high statistics MC signal event sample generated according to a double-differential 
cross section d 2 a gener i C / (dx dQ 2 ) obtained from (1) by replacing the Breit-Wigner LQ propaga- 
tor term with a constant. This unique MC sample is used to calculate the efficiency to select a 
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LQ of a given type, mass m L Q, coupling \ eq and /?lfv by attributing to each event a weight: 



w(type, m LQ , A eg , /3 LF v, z, <? 2 ) 




(type, m LQ , A eg , /3 LFV , sc, Q 2 ) 



(3) 




where Q 2 = sxy refers to the generated negative momentum transfer squared and x is the 
Bj0rken scaling variable known at the generator level. This procedure provides an exact pre- 
diction over the full range of LQ production parameters and avoids approaches like the narrow 
width approximation or the high mass (contact interaction) approximation. 

The LQ kinematics are reconstructed using the double angle method [16]. The direction of 
the detected lepton and jet are used to reconstruct the Bj0rken scaling variable x and therefore 
the LQ mass m^Q = y/xs. 

The contributions from Standard Model (SM) background processes which may mimic the 
signal include neutral current (NC) and charged current (CC) deep-inelastic scattering (DIS), 
photoproduction, lepton pair production and real W boson production. These processes are 
briefly described below: 

• NC DIS (ep -> eX) 

NC DIS processes contribute to the selected event sample if the scattered electron is 
attributed to the tau electronic decay or if it is misidentified as a narrow jet corresponding 
to a tau decay to hadrons. The NC DIS background is modelled using the event generator 
RAPGAP [17]. The proton PDFs are parametrised using CTEQ5L [15] and hadronisation 
is performed using JETSET [18] parton showers and the Lund string fragmentation. 

• CC DIS (ep -> uX) 

Lepton flavour violating processes usually exhibit an imbalance in the measured calori- 
metric transverse momentum due to either the presence of a minimally ionising muon 
in fiX final states or the escaping neutrino(s) from tau decays in rX events. This im- 
balance is exploited in the event selection. The CC DIS process leads to events with 
genuine missing transverse momentum and therefore contributes to the selected sample 
if hadrons or photons from the final state are misidentified as muons or if tau decays are 
falsely reconstructed. The CC DIS contribution is modelled using the DJANGO event 
generator [19]. 

• Photoproduction (77? — > X) 

Events from photoproduction processes may contribute to the final selection if a hadron 
is wrongly identified as a muon or if a narrow hadronic jet fakes the tau signature. This 
contribution is calculated using the event generator PYTHIA [20]. CTEQ5L [15] serves 
as the proton PDF parametrisation and the photonic parton distribution parametrisation 
GRV-LO [21] is used. As PYTHIA only contains leading order 2 — > 2 processes, the 
multi-jet production cross section is underestimated [22]. Therefore, the prediction is 
scaled up by a factor 1.2 in this analysis, in agreement with previous analyses of jets in 
photoproduction [22] . 
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• Lepton-pair production (ep — > e£ + £~X) 

Lepton-pair production events contribute to the background because they may lead to 
high momentum leptons in the final state. In particular, inelastic di-muon events with one 
unidentified muon may fake the fxX LFV signature. The background samples include ee, 
fx fx and rr production generated with the event generator GRAPE [23]. 

• W production (ep — > eWX) 

Real W boson production leads to final states with isolated high P T leptons and missing 
transverse momentum. The simulated W production samples are created with the event 
generator EPVEC [24] and include leptonic (eu e , fxu^, tu t ) and hadronic W decays. 

All signal and SM samples are passed through a detailed simulation of the HI detector 
response based on the GEANT program [25] and the same reconstruction and analysis algo- 
rithms as used for the data. 

4 High Pt Muon Signatures 

Leptoquarks with couplings to the first and the second lepton generation can be produced in ep 
collisions and may decay to a muon and a quark. The signature is an isolated high Pt muon 
back-to-back to the hadronic system in the transverse plane. In general, a muon deposits a very 
small fraction of its energy in the LAr calorimeter. The signal is therefore expected to exhibit 
large P£ al °, which is the net transverse momentum reconstructed from all clusters recorded in 
the LAr calorimeter alone. 

The event preselection requires at least one muon with a transverse momentum above 10 GeV 
in the polar angular range 10° to 140° and at least one jet. The muon is required to be isolated. 
The angular distance, D = \J (Arj) 2 + (A0) 2 , of the muon to the nearest track and to the near- 
est jet is required to be greater than 0.5 and 1.0, respectively. Only events with P£ al ° greater 
than 12 GeV are selected. In order to further exploit the event topology in the transverse plane, 
the cut V ap /Vp < 0.3 is employed, where V ap /V p is defined as the ratio of the anti-parallel to 
parallel projections of all energy deposits in the calorimeter with respect to the direction of 
P£ al ° [26]. 

Figure 2 displays the distributions of the transverse momentum of the muon, its polar angle 
Op, P T al ° and the acoplanarity A0 M _ x between the muon and the hadronic final state X after 
the muon preselection. The data passing the preselection are well described by the SM predic- 
tion. The signal corresponding to a scalar LQ with m LQ = 200 GeV is also shown. It displays 
muons with large P^ produced predominantly in the forward direction (low 6^) in events with 
significant P£ al ° and back-to-back topology A0 M _x — 180°. 

In the final LFV selection step, the NC DIS background is further suppressed by reject- 
ing events with identified electrons, and by accepting only events with an imbalance of the 
calorimeter deposits, P£, al ° > 25 GeV, and with a back-to-back topology, A0 M _x > 170°. The 
latter selection criterion is only applied for events for which the hadronic final state is well 
contained in the detector, with the reconstructed polar angle 7° < 9\ < 140°. 

The selection efficiency ranges from 40% to 60% depending on the LQ mass and type (see 
table 2). 
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5 High P T Tau Signatures 



Leptoquarks with couplings to the first and the third lepton generation can be produced in ep 
collisions and may decay to a tau and a quark. Tau leptons are identified using the electronic, 
muonic and hadronic decays of the tau. 

5.1 Electronic tau decays 

The final state resulting from the electronic tau decay, r — > ev e v T , leads to an event topology 
that is very similar to that of high Q 2 NC DIS events. The preselection follows that presented 
in [27]. A reconstructed jet with a minimal transverse momentum of Pj, > 25 GeV back- 
to-back in the transverse plane to an electron with Pj, > 10 GeV is required. The kinematic 
domain is restricted to Q 2 > 1000 GeV 2 and y > 0.1. Figure 3(a) shows the distribution of 
P™ 1SS after this preselection, where P™ 1SS is defined as the total missing transverse momentum 
reconstructed from all observed particles. 

In the final selection a large missing transverse momentum P™ 1SS > 20 GeV is required in 
order to account for the expected missing momentum carried by the neutrinos produced in the 
tau decay. These neutrinos are boosted along the electron direction, implying an imbalance 
between the transverse momenta of the electron Pj, and the hadronic final state P* . Hence, 
the restriction Pf/P^ < 0.8 further reduces NC DIS background. In addition, the azimuthal 
distance between the missing transverse momentum and the electron must not exceed 20°. The 
remaining NC DIS background, due to mismeasured electron energies leading to missing en- 
ergy near the electron, is reduced by the requirement P^~ clu / P^~ trk > 0.7, where P^~ cln is 
measured from the electromagnetic cluster and P|T trk from tracking information. Any events 
with additional isolated muons are excluded from the electronic tau decay channel. The final 
selection in the electronic tau decay channel yields an efficiency normalised to all tau decays 
of 3% to 10%, which is limited by the branching fraction BR(r — > ev e v T ) = 17.8% [28] and 
dependent on the assumed LQ mass and type (see table 2). 

5.2 Muonic tau decays 

Muonic tau decays r — »■ nv^v T result in similar final states as the high Pr muon signatures 
described in section 4. The same selection cuts described therein are applied here. To ac- 
count for possible effects due to different muon kinematics resulting from a tau decay, the 
selection efficiency was studied in detail with a LFV MC signal sample with a rX final state 
and a subsequent muonic tau decay. The selection efficiency ranges between 4% and 8%, 
which is dependent on the LQ mass and type, normalised to all tau decays and limited by 
BR(t -> fiu^Vr) = 17.4% [28] (see table 2). 

5.3 Hadronic tau decays 

The hadronic decays of the high Pr tau lead to a typical signature of a high Pr "pencil-like" jet. 
The signal topology is a di-jet event with no leptons. The tau-jet is characterised by a narrow 
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energy deposit in the calorimeter and low track multiplicity with predominantly one or three 
tracks in the identification cone of the jet. The neutrino from the tau decay are boosted along 
the direction of the hadrons. The missing transverse momentum in the event is aligned with the 
tau-jet. 

Tau-jet candidates are defined as jets with exactly one or three tracks in the tau-jet cone 
with an opening angle that varies between 5° and 30° with decreasing jet momentum. The 
tracks are required not to be associated with identified electrons or muons and the scalar sum 
of their transverse momenta is required to be larger than 2 GeV. The fine granularity of the LAr 
calorimeter is used to match extrapolated tracks with energy deposits in the calorimeter and 
to separate additional neutral particles associated to the tau candidate from unmatched energy 
deposits in the tau-jet cone. The sum of the four-vectors of the tracks and of the neutral particles 
defines the tau-jet candidate four-vector. 

In the preselection step at least two jets with a transverse momentum Pj, otl > 20 GeV and 
P^ et2 > 15 GeV reconstructed in the polar angle range 7° < #j ets < 145° are required. One jet 
must fulfil the criteria of a tau-jet candidate with 9 T j Ct > 20°. In addition, the calorimetric 
shower shape and tracking signature are exploited to validate the tau-jet candidates. The fol- 
lowing estimators are used to separate a tau-jet from quark or gluon induced jets: the number 
of all tracks associated to the tau-jet candidate, the distance in rj — (ft between tracks and calori- 
metric clusters, the number of calorimeter cells of the tau-jet n ce n s , the radial extension of the 
calorimetric deposits (r) = ^"f 6 " 8 Ei r i/ Yli-Ei, the standard deviation a(r) = (r 2 ) — (r) 2 
and the invariant tau-jet mass reconstructed from calorimeter cells. A neural net algorithm is 
employed and trained using the six estimator variables, as explained in [29]. The neural net 
yields a discriminator variable V NN in the range < Pnn < 1 with values close to for quark 
or gluon induced jets and close to 1 for hadronic tau decays. The distribution of the discriminant 
Vnn after the preselection is depicted in figure 3(b). The distributions of P™ 1SS and A0 miss „ r j Ct 
after requiring V^n > 0.8 are shown in figures 3(c),(d). This requirement yields a signal effi- 
ciency of 80% and a quark or gluon induced jet rejection of 95%. After all preselection criteria 
16 (112) events are selected in e~p (e + p) data sample for 22.0 ± 1.0(stat.) (121.1 ± 5.3(stat.)) 
expected from the SM. 

The final selection step in the hadronic tau decay channel makes use of the characteristic 
large missing transverse momentum carried by the tau neutrino which is expected to be in the 
direction of the tau-jet. The difference in <p between the missing transverse momentum vector 
and the tau-jet, A0 miss _ T j et , is required to be below 20°. A minimal value of P™ 1SS > 12 GeV 
is chosen for an accurate determination of the direction. In addition P£, al ° > 12 GeV is re- 
quired. The final signal selection efficiency in the hadronic tau decay channel varies between 
3% and 13%, normalised to all tau decays and limited by the branching fraction BR(t~ — > 
v T + hadrons) = 64.8% [28] (see table 2). 



6 Systematic Uncertainties 

The following experimental systematic uncertainties are considered: 
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• The energy of electrons is measured with a systematic uncertainty in the range from 0.7% 
to 3% depending on the polar angle. The uncertainty of the electron direction is estimated 
to be less than 3 mrad in 9 and 1 mrad in 0. 

• The scale uncertainty on the transverse momentum of high Pt muons amounts to 5%. 
The uncertainty on the reconstruction of the muon direction is 3 mrad in 9 and 1 mrad in 
0. 

• For the hadronic final state, an energy scale uncertainty of 2% and a direction uncertainty 
of 20 mrad are assumed. 

• The luminosity of the analysed datasets is known to 1.5%. 

The effects of these systematic uncertainties on the signal and the expected SM background 
are evaluated by shifting the relevant quantities in the MC simulation by their uncertainty and 
adding all resulting variations in quadrature. 

Systematic errors accounting for normalisation uncertainties on the expected background 
determined from the individual MC event generators are estimated to be 10% for NC DIS and 
Lepton-pair production, 15% for W production and 30% for photoproduction and CC DIS. 
The relatively large error of 30% on photoproduction and CC DIS is due to uncertainties on 
higher-order corrections. The errors associated to the background normalisation are added in 
quadrature to the experimental error to calculate the total error of the SM prediction. 

The main theoretical uncertainty on the signal cross section originates from the parton den- 
sities. This uncertainty is estimated as described in [13]. It is found to be 5% for LQs coupling 
to up-type quarks and varies between 7% at low masses and 30% at masses around 290 GeV 
for LQs coupling to down-type quarks. The correlation between different channels is taken into 
account for the statistical interpretation and limit calculation [30]. A detailed description of the 
analysis can be found in [31]. 

7 Results 

No candidate is found in the final data sample of the muon channel. The expected number of SM 
background events is 1.03 ± 0.32 in the e + p set and 0.18 ± 0.06 in the e~p sample. The largest 
contribution to this background comes from muon-pair production and the muonic decays of 
W bosons. These results apply equally to the muonic tau decay channel. 

In the electronic tau decay channel no data event is found compared to a SM expectation of 
0.28 ± 0.19 events in the e~p sample and 1.24 ± 0.55 events in the e + p data. NC DIS events 
with a mismeasured electron energy are the largest background contribution. 

No e~p data event passes the final selection criteria in the hadronic tau decay channel. The 
expected SM background amounts to 0.29 ± 0.06. One event is selected in the e + p data for an 
expected SM prediction of 2.63 ± 0.57, dominated by NC DIS and photoproduction processes. 

The results of the final selection in all channels are summarised in table 2. Typical signal 
selection efficiencies for some LQ types with a mass of 150 GeV and 500 GeV are also given. 
The observation is in agreement with the SM prediction and no evidence for LFV is found by 
the present analysis. Limits on the model parameters presented in section 3 are calculated as 
described in the following section. 
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8 Limits 



The results of the search are interpreted in terms of exclusion limits on the mass and the coupling 
of LQs that may mediate LFV. The LQ production mechanism at HERA involves non-zero 
coupling to the first generation fermions X eq > 0. The LFC decays LQ— >eg or LQ^u e q are 
therefore possible. In order to cover the full LQ decay width and to generalise the results of LFV 
searches in ep collisions to an arbitrary weight between the LFC and LFV decay channels, the 
searches for LFC decays presented in [13] are combined with each of the LFV search channels 
fj,X or tX of the present analysis. It is assumed that only one of the couplings X m and X Tq 
is non-zero and therefore the fxX and tX channels enter the limits calculation separately. A 
modified frequentist method with a likelihood ratio as the test statistic is used to combine the 
individual data sets and search channels [32]. 

In first generation LQ signals are searched for in about 400 bins in the m L q — y plane and 
the observed data is in agreement with the irreducible SM NC and CC background [13]. For 
the LFV channel [iq (rq), the couplings X eq and X m (X Tq ) and the LQ mass determine the 
total production cross section, which is compared to the selected data from the LFV search 
channel and the first generation results. A combined test statistic is built and used to set limits 
as a function of X eq , A Aig (A rg ) and m L Q. This procedure implicitly includes in the analysis the 
decays to a neutrino of any flavour and a quark. 

Figure 4 shows limits before and after combination with the search for first generation LQs 
for the LQ types and V L up to LQ masses of 320 GeV assuming X eq = X m and X eq = X Tq , 
i.e. Plfv = 0.5, in the resonance production region. The comparison for these types exemplifies 
that the limits on those LQs which can decay to a neutrino-quark pair, namely Sq , S'f , Vq and 
Vi, benefit most from the combination with the search for first generation LQs which covers 
decays to a neutrino-quark pair. In the high mass regime m L q \fs (contact interaction region) 
the obtained limits are similar to those deduced without the combination. The fluctuations in 
the combined limits are due to the observed data events in the search for first generation LQs. 
In the mass range from 250 GeV to 300 GeV both the combined limits on X m and X Tq are for all 
LQ types up to a factor 2 more stringent than without combination. Table 3 shows the 95% CL 
combined lower limits on the LQ mass for all LQ types assuming a coupling of electromagnetic 
strength X eq = X M (X Tq ) =0.3. 

Allowing for an arbitrary decay rate between the LFC and LFV decay channels, /3lfv> the 
excluded regions for two LQ types and four mass values in the X Ml — X eqi (a,b) and X rqi — X eqi 
(c,d) planes are presented in figure 5. For very low values of /?lfv (A eg 3> X M (X rq )), the limits 
on X eq reproduce the bounds published in [13], as expected, since the LFC channel dominates 
the LQ width. For /^lfv^O.S (A w (A T(? ) ^> X eq ) the present analysis extends significantly the 
published limits on X eq to lower values. The limit without combination in the contact interac- 
tion region (where the cross section is proportional to X eqi X^ T ) q Jm^^) forms a cross-diagonal 
straight line following different values of Plfv- The combination in the contact interaction re- 
gion, e.g. mLQ = 350 GeV, barely strengthens the limit as the virtual effects of the high mass 
LQ contact interaction at low values of \/I are marginal compared to the irreducible NC and CC 
DIS background. Fluctuations of the data may even result in a less stringent combined limit. 

Figures 6 and 7 display the 95% CL upper limits on the coupling X m and X Tq of all 14 LQ 
types to a muon-quark pair and a tau-quark pair, respectively, as a function of the LQ mass 
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leading to LFV in ep collisions, assuming X eq = A M (A Tg ) The limit curves referring to the LQ 
types Sq and SK 2 are identical to profiles of the corresponding excluded regions following 
the value /?lfv = 0.5 in figure 5. The limits are most stringent at low LQ masses with values 
C(10 -3 ) around ttilq = 100 GeV. Corresponding to the steeply falling parton density function 
for high values of x, the LQ production cross section decreases rapidly and exclusion limits 
are less stringent towards higher LQ masses. For LQ mass values near the kinematical limit of 
319 GeV, the limit corresponding to a resonantly produced LQ turns smoothly into a limit on 
the virtual effects of both an off-shell s-channel LQ process and a u-channel LQ exchange. At 
masses m L Q > \fs the two processes contract to an effective four-fermion interaction, where the 
cross section is proportional to (\^ q \ eq / m^^) 2 . This feature is visible in the constant increase 
of the exclusion limit for masses above the ep centre-of-mass energy of yfs = 319 GeV. Due to 
initial state QED radiation and very low parton densities for masses near s/s the "kink" of the 
transition region is shifted to somewhat smaller masses of around 290 — 300 GeV. 

It is noticeable that the limits on vector LQs are more stringent compared to those on the 
scalars, due to the considerably higher cross section and the slightly higher acceptance. In each 
plot those LQ types that have couplings to both u and d quarks exhibit the best limit. The limits 
corresponding to LQs coupling to a u quark are more stringent than those corresponding to 
LQs coupling to the d quark only, as expected from the larger u quark density in the proton. 
The LQs Sq and S R (V l and V^) differ only by the decay into a neutrino and a quark of the 
lefthanded LQ. As this decay channel is not covered in the LFV decay channels, the left-handed 
LQ cannot be as strictly excluded as the right-handed one. This argument applies to the resonant 
production where the analysis is only sensitive to the partial width of the LQ. In the high mass 
region the limits for Sq and Sq (V l and V R ) are similar, as the four-fermion interaction is 
independent of the decay width. 

The limits on \^( T ) q = \ eq derived from the virtual effects of a 500 GeV LQ are transformed 
into a limit on the value A M ( T )^ A e9i /mLQ and shown in tables 4 and 5 for F = LQs and in 
tables 6 and 7 for F = 2 LQs. For each LQ type the limit is calculated for the hypothesis of 
a process with only the quarks of flavours i and j involved. With respect to quark flavours, 
the selection criteria described in sections 4 and 5 are inclusive since no flavour tagging of the 
hadronic jet is used. 

These results may be compared with constraints from low energy experiments, based on 
the non-observation of LFV in muon scattering and rare decays of mesons and leptons [28]. 
The interpretation in terms of leptoquark exchange and limits on A ft ( r ) (?j -Ae^/m^g [33] are also 
shown in tables 4, 5, 6 and 7. Bounds of similar magnitude are observed for processes involving 
e — > r transitions and charm or bottom quarks. In these cases the limits obtained in the present 
analysis are often superior to those from low energy experiments. 

The results on LFV in LQ production are directly comparable with those from the ZEUS 
experiment [34]. Similar limits are obtained. At hadron colliders LQs are mainly produced in 
pairs independently of the coupling, and therefore searches cannot constrain LFV couplings. 
Lower mass limits on the second and third generation leptoquarks extend up to 250 GeV and 
150 GeV, respectively, depending on the type and the assumed decay branching ratios [35]. 
Similarly, second and third generation leptoquarks are pair produced in e + e~ annihilation where 
typical lower mass bounds reach values of 100 GeV [36]. 
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9 Conclusion 



A search for lepton flavour violation processes induced by leptoquarks in ep collisions at a 
centre-of-mass energy of 319 GeV with the HI experiment at HERA is presented. No signal for 
the LFV processes ep — > fxX or ep — > rX is found. Constraints on LFV LQ couplings are set 
combining the LFV search with the search for first generation LQs. The limits are a factor of 2 
to 4 more stringent and extend beyond the domain in LQ mass excluded by previous searches 
performed by the HI experiment [7]. Exclusion limits on several scenarios of LFV transitions 
of the kind eq,i — > rqj are more stringent than limits from searches for certain rare meson or tau 
decays. Assuming a coupling of electromagnetic strength, leptoquarks mediating lepton flavour 
violating processes e — > fi and e — > r can be ruled out up to masses of 459 GeV and 379 GeV, 
respectively. 
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Table 1 : The 14 leptoquark (LQ) types of the Buchmuller-Riickl-Wyler classification [12] in the 
Aachen notation [14]. The LQ subscripts refer to the weak isospin and the superscripts refer to 
the lepton chirality. Columns 2-4 display the spin J, fermion number F and electrical charge 
Q. The dominant resonant production process in ep scattering and the corresponding coupling 
is shown in columns 5 and 6 respectively. Leptoquarks which couple to a left-handed lepton 
doublet and can decay into a neutrino-quark pair, have a charged lepton decay branching ratio 

fA=rv(r^ + r^) = i/2. 
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Table 2: Summary of the selection results of the search for the LFV processes ep — > [iX and 
ep — > tX. The results of the individual tau decay channels are also shown. The errors on the 
SM MC expectation include statistical and systematic errors added in quadrature. Examples 
of signal selection efficiencies for leptoquarks of the types S$, Vy 2 , and Sy 2 with masses 
mLQ of 150 GeV and 500 GeV are also shown. For the tau decay channels the efficiencies are 
normalised to the sum of all tau decays. 
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F = 2 




qR 


4* 




V l/2 


>/ l/2 


V L 

V l/2 


eq — > /xg 


294 


294 


278 


306 


299 


374 


336 


eg — > rg 


293 


294 


276 


295 


282 


302 


297 



Table 3: Lower exclusion limits at 95% CL on leptoquark masses m L Q assuming 

\t<? = ^eq — 0.3 or \ Tq = X eq = 0.3. 
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ep - 






HI 




F 


= 




Upper exclusion limits on A e9l A 


w >lq (TeV- 2 ) 




for lepton flavour violating leptoquarks at 95% CL 




^1/2 

e~ u 


oR 
°l/2 

e~ {u + d) 
e+(u + d) 


^1/2 

e~ d 
e+d 


v 

e~ d 
e+d 


K 

e~ d 
e+d 


Vn R 

v 

e~ u 
e~^~u 


y 1 

e~(V2u + d) 
e+iVZu, + d) 


1 1 


[AN -> eN 
5.2 x 10 -5 

1.4 


fiN -> eJV 
2.6 X 10~ 5 
1.1 


fiN -> eiV 
5.2 x 10~ 5 

2.0 


fiN eiV 
2.6 X 10~ 5 

1.3 


fiN -> eN 
2.6 X 10~ 5 

1.3 


fiN — eN 
2.6 X 10~ 5 

0.9 


fiN -» eN 
0.8 X 10~ 5 

0.4 


1 2 


D — * fie 
2.4 


ff — ► ^ie 
2 x 10~ 5 


K — > fie 
2 x 10~ 5 


K - -> jj,e 
1 X 10 -5 


K -> fie 
1 X 10~ 5 


D — ► fie 
1.2 


K — > fie 
1 x 10~ 5 




1.4 


1.2 


2.0 


1.5 


1.5 


1.1 


0.5 


1 3 


* 


B — > fie 
0.4 

2.1 


B — ► fie 
0.4 

2.1 


B — v 
0.2 

1.6 


B — * /j,e 
0.2 

1.6 


* 


B — y fie 
0.2 

1.6 


2 1 


D — > fie 
2.4 

4.2 


K — > ^ie 
2 x 10~ 5 

2.9 


K — > fie 
2 x 10~ 5 

4.1 


/C — ► /^e 
1 X 10~ 5 

1.7 


K — ► fie 
1 x 10~ 5 

1.7 


D — ► fie 
1.2 

1.5 


K — > fie 
1 X 10~ 5 

0.7 
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MiV eiV 
9.2 x 10~ 4 

6.0 


fiN -> eJV 
1.3 X 10~ 3 

3.7 


fiN -> eJV 
3 x 10~ 3 

4.8 


fiN eiV 
1.5 X 10~ 3 

2.5 


fiN -» eJV 
1.5 X 10~ 3 

2.5 


fiN -> eJV 
4.6 X 10 -4 

3.1 


fiN -> eAT 
2.7 X 10~ 4 

1.3 


23 


* 


0.3 

5.2 


B — * fieK 
0.3 

5.2 


B —v /ieif 
0.15 

3.5 


B -> peJC 
0.15 

3.5 


* 


B -» fieA' 
0.15 

3.5 


3 1 


* 


B — ► fie 
0.4 

5.3 


B — * fie 
0.4 

5.3 


0.12 

1.8 


_B — ► fie 
0.2 

1.8 


* 


0.12 

1.8 


32 


* 


B -> fieK 
0.3 

7.0 


B -t fieK 
0.3 

7.0 


S —y fieK 
0.15 

2.8 


B -» fieK 
0.15 

2.8 


* 


B -> fieK 
0.15 

2.8 


33 


* 


/jiV -t eN 
1.3 X 10 -3 

8.3 


fiN -> eN 
3 x 10~ 3 

8.3 


/jJV -> eJV 
1.5 x 10~ 3 

4.3 


fiAT -» eJV 
1.5 x 10 -3 

4.3 


* 


fiN -» eAT 
2.7 x 10" 4 

4.3 



Table 4: Limits at 95% CL on A^A^ /m^n for F = leptoquarks (bold). Combinations of i 
and j shown in the first column denote the quark generation coupling to the electron and muon 
respectively. In each cell the first two rows show the process providing the most stringent limit 
from low energy experiments. The cases marked with '*' refer to scenarios involving a top 
quark, not considered in the present analysis. Highlighted HI limits are more stringent than 
those from the corresponding low energy experiment. 
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ep - 






HI 




F 


= 




Upper exclusion limits on A e?4 A 


rqj /mt Q (TeV~ 2 ) 


for lepton flavour violating leptoquarks at 95% CL 




oL 
u l/2 

e~ u 
e+ u 


qR 

°l/2 

e~ (u + d) 
e+(u + d) 


oL 
u l/2 

e~ d 
e+d 


v 

e~d 
e+d 


y 

e~ d 
e+d 


K 

e — li 
e + u 


vf 

v 1 

e~ (VM + 3) 
e+fv^u + d) 


1 1 


r — > 7re 
0.4 

2.1 


r — > 7re 
0.2 

1.8 


T — ► rre 
0.4 

3.1 


T — > rre 
0.2 

2.1 


r — ' rre 
0.2 

2.1 


r — > rre 
0.2 

1.5 


r — * 7re 
0.06 

0.7 


1 2 




r — * Ke 
6.3 


K — > ■kuu 
5.8 X 10~ 4 


T ^ Ke 
3.2 


r — * Ke 
3.2 




K — > nfis 
1.5 X 10~ 4 




2.2 


1.8 


3.2 


2.5 


2.5 


1.7 


0.8 


1 3 


* 


B -> re 
0.3 

3.2 


0.3 

3.2 


B — > re 
0.13 

2.7 


B -» re 
0.13 

2.7 


* 


B -> re 
0.13 

2.7 


2 1 






K — ► TYUU 


r-tJTe 


r -» if e 




K — * -kvv 




6.3 


5.8 X 10 -4 


3.2 


3.2 




1.5 X 10 -4 




6.7 


4.8 


6.9 


2.8 


2.9 


2.3 


1.1 
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t — * 3e 


r — + 3e 


r — > 3e 


r — • 3e 


r — + 3e 


r — ► 3e 


r — > 3e 


5.0 


8.0 


17.0 


9.0 


9.0 


3.0 


1.6 




10.9 


6.7 


8.6 


4.5 


4.5 


5.5 


2.4 


23 


* 


B -> reX 
14.0 


B -> reX 
14.0 


B -> rel 
7.2 


B -► reX 


* 


B -► reX 
7.2 






9.3 


9.3 


6.3 


6.3 




6.3 


3 1 


* 


B —> TE 

0.3 

9.1 


B -> re 
0.3 

9.1 


0.12 

3.0 


B re 
0.13 

3.0 


* 


0.12 

3.0 


32 


* 


B -> reX 
14.0 


B -i reX 
14.0 


B -► rel 
7.2 


B -> reX 
7.2 


* 


B -► rel 
7.2 






12.6 


12.6 


4.9 


4.9 




4.9 


33 


* 


r — * 3e 
8.0 


t — > 3e 
17.0 


r — ' 3e 
9.0 


r — ' 3e 
9.0 


* 


r — ► 3e 
1.6 






15.2 


15.2 


8.1 


8.1 




8.1 



Table 5: Limits at 95% CL on \ eqi \ Tq J for F = leptoquarks (bold). Combinations of 
i and j shown in the first column denote the quark generation coupling to the electron and tau 
respectively. In each cell the first two rows show the process providing the most stringent limit 
from low energy experiments. The cases marked with '*' refer to scenarios involving a top 
quark, not considered in the present analysis. Highlighted HI limits are more stringent than 
those from the corresponding low energy experiment. 
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ep -> fiX HI F = 2 


Upper exclusion limits on Ae^A^/m^g (TeV 2 ) 
for lepton flavour violating leptoquarks at 95% CL 


QiQj 


^0 

e — u, 




e — u 




e~ (it + d) 
e+(u + d) 


e~ (u + v'Sd) 
e+(it + V2d) 


K 1/2 

e — d 
e+d 


v 1/2 

e~ (u + d) 
e+(u + d) 


v 1/2 

e — iz 


1 1 


[AN -> eN 
5.2 x 10~ 5 

2.0 


MiV -* eiV 
5.2 x 10 -5 

2.0 


/jJV -> eJV 
5.2 x 10~ 5 

2.6 


liN -> eiV 
1.7 x 10~ 5 

1.0 


liN -> eN 
2.6 X 10~ 5 
1.1 


liN -> eJV 
1.3 X 10~ 5 

0.6 


f*iV -» eAT 
2.6 X 10~ 5 

0.8 


1 2 


K — ► 1YVU 

1 x 10 -3 

2.6 


£> -» fie 
2.4 

2.6 


K — > fie 

2 x 10 -5 

3.2 


A — ► fie 
1 X 10 -5 

1.4 


K — v fie 
1 x 10 -5 

2.0 


K — > fie 
1 x 10 -5 

1.4 


Z) — > /j,e 
1.2 

1.9 


1 3 


* 


* 


S — ► fie 
0.4 

3.3 


0.24 

1.6 


B — ► fie 
0.2 

2.5 


B — ► fie 
0.2 

2.5 


* 


2 1 


K — * tyuu 
1 X 10 _ 3 

2.6 


D — > fie 
2.4 

2.6 


A' — > fie 
2 x 10~ 5 

3.3 


/C — > fie 
1 X 10~ 5 

1.4 


K — ► fie 
1 x 10~ 5 

1.1 


K — > fie 
1 x 10~ 5 

0.7 


D — > /j,e 
1.2 

0.8 
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MiV eN 
9.2 x 10~ 4 

6.0 


/jiV -> eN 
9.2 x 10~ 3 

6.0 


fiN -> eJV 
3 x 10~ 3 

4.8 


fiN -> eAT 
2.5 X 10 -3 

2.2 


fiN -» eiV 
1.5 X 10~ 3 

2.5 


fiN -> eiV 
6.7 x 10~ 4 

1.9 


fiN -t eiV 
4.6 X 10~ 4 

3.1 


23 


* 


* 


B -> peA" 
0.3 

5.2 


B — » /ieA' 
0.15 

2.6 


B -> fieK 
0.15 

3.5 


B -> fieK 
0.15 

3.5 


* 


3 1 


* 


* 


S — * fie 
0.4 

3.7 


B —• lie 
0.4 

1.9 


B — ► fie 
0.2 

1.2 


B — * fie 
0.2 

1.2 


* 


32 


* 


* 


B -t fieK 
0.3 

7.0 


S — <■ jleK 
0.15 

3.5 


B -► fieK 
0.15 

2.8 


B — fieK 
0.15 

2.8 


* 


33 


* 


* 


fiN -> eJV 
3 X 10~ 3 

8.3 


liN -> eiV 
2.5 X 10 -3 

4.3 


fiAT — * eJV 
1.5 X 10 -3 

4.3 


f,N -> eN 
6.7 x 10~ 4 

4.3 


* 



Table 6: Limits at 95% CL on A^A^ /m^n for F = 2 leptoquarks (bold). Combinations of i 
and j shown in the first column denote the quark generation coupling to the electron and muon 
respectively. In each cell the first two rows show the process providing the most stringent limit 
from low energy experiments. The cases marked with '*' refer to scenarios involving a top 
quark, not considered in the present analysis. Highlighted HI limits are more stringent than 
those from the corresponding low energy experiment. 
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ep -> tX HI F = 2 



Upper exclusion limits on X eqi X Tqj /m1q (TeV 2 ) 
for lepton flavour violating leptoquarks at 95% CL 





e~ u 
e~^~u 


e — u 
e 4 " u 
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e~ (u + ^/2d) 
e+(u + \/2d) 


1/2 

e — d 
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e~ (it + d) 
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r — ► 7re 
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3.0 


0.4 

3.0 


0.4 
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0.1 

1.7 


0.2 

1.8 


0.1 

1.0 


0.2 

1.2 


1 2 


K — * TtVV 




r -> Ke 


— > TZVV 


K — > nut/ 
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5.8 X 10~ 4 




6.3 


2.9 x 10~ 4 


2.9 X 10~ 4 


3.2 






4.0 


4.0 


5.0 


2.1 


3.5 


2.3 


3.1 


1 3 


* 


* 


fl -t re 
0.3 

5.3 


0.12 

2.7 


B -» re 
0.13 

4.2 


B -> re 
0.13 

4.2 


* 


2 1 


K — * -kvv 






K — > 1YVU 


K — ► irviy 


r — > i-v e 




5.8 X 10~ 4 




6.3 


2.9 x 10 -4 


2.9 X 10~ 4 


3.2 






4.2 


4.2 


5.5 


2.3 


1.8 


1.1 


1.2 


2 2 


r — * 3e 


r — » 3e 


r — * 3e 


t — * 3e 


r — ► 3e 


r — ► 3e 


r — )■ 3e 


5.0 


5.0 


17.0 


14.0 


9.0 


4.0 


3.0 




10.8 


10.9 


8.6 


3.9 


4.5 


3.5 


5.5 
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* 


B -> feX 
14.0 


B -> reX 
7.2 


B reX 


B -> reX 
7.2 


* 








9.3 


4.7 


6.3 


6.3 




3 1 


* 


* 


B -» re 
0.3 

6.3 


B — > re 
0.13 

3.1 


B -> re 
0.13 

1.9 


B -► re 
0.13 

1.9 


* 


32 


* 


* 


B -» feX 
14.0 


B -> reX 
7.2 


B -► reX 
7.2 


B -> feX 
7.2 


* 








12.6 


6.4 


4.9 


4.9 




3 3 


* 


* 


r — ► 3e 
17.0 


T -> 3e 
14.0 


t — > 3e 
9.0 


r — )■ 3e 
4.0 


* 








15.2 


7.8 


8.1 


8.1 





Table 7: Limits at 95% CL on \ eqi \ Tq J for F = 2 leptoquarks (bold). Combinations of 
i and j shown in the first column denote the quark generation coupling to the electron and tau 
respectively. In each cell the first two rows show the process providing the most stringent limit 
from low energy experiments. The cases marked with '*' refer to scenarios involving a top 
quark, not considered in the present analysis. Highlighted HI limits are more stringent than 
those from the corresponding low energy experiment. 
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Figure 2: Control distributions of the preselected fiX sample: (a) muon transverse momentum, 
(b) muon polar angle, (c) transverse momentum as measured from the calorimeter deposits and 
(d) acoplanarity between the muon and the hadronic final state X. Data (points) from e + p 
collisions are compared to the SM expectation (histogram). The LFV signal MC sample of a 
leptoquark Sy 2 with m L Q = 200 GeV and X eq = X m = 0.3 is shown hatched with arbitrary 
normalisation in each plot. 
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HI 
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Figure 3: Distributions of the preselected rX sample: (a) missing transverse momentum in the 
electronic tau decay channel and (b) neural net tau-jet discriminant after the preselection in the 
hadronic tau decay channel. The restricted sample obtained after the additional cut X> NN > 0.8 
in the hadronic channel: (c) missing transverse momentum and (d) acoplanarity between the 
tau-jet and the missing transverse momentum. The LFV signal MC sample of a leptoquark Sy 2 



with m LQ 
each plot. 



200 GeV and A 



eq 
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0.3 is shown hatched with arbitrary normalisation in 
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Figure 4: Comparison of limits at 95% CL on the coupling constants \g q under the assumption 
Afo = X eq as a function of the leptoquark mass m L Q for: (a) Sfr on A MU = X eu , (b) V L on X,,a = \„ 



(c) Sq on X TU = \ eu , and (d) Vq on \ T d = A e d. The areas above the dashed lines represent the 
exclusion regions using only the lepton flavour violating leptoquark decay channels ep — > /xX 
and ep — > rX, respectively. The limits after combination with the results of the search for first 
generation leptoquarks are shown as solid lines. 
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HI Search for lepton flavour violation 




Figure 5: Excluded regions at 95% CL (filled) on Xc qi as a function of X eqi for four differ- 
ent leptoquark masses. The branching ratio /?lfv = ^e qi /(^j qi + ^l qi ) is not fixed. Diagonal 
dashed lines represent iso-curves for fixed values of /3lfv- The bounds deduced without the 
combination with first generation leptoquarks are shown as black curves corresponding to the 
different mass assumptions. 
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HI Search for lepton flavour violation 
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Figure 6: Limits on the coupling constants X Ml = \ eqi as a function of the leptoquark mass m L Q 
for (a), (b) F — and (c), (d) F = 2 scalar and vector leptoquarks. Regions above the lines are 
excluded at 95% CL. The notation q 1 illustrates that only processes involving first generation 
quarks are considered. 
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HI Search for lepton flavour violation 
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Figure 7: Limits on the coupling constants X rqi = \ eqi as a function of the leptoquark mass m L Q 
for (a), (b) F = and (c), (d) F = 2 scalar and vector leptoquarks. Regions above the lines are 
excluded at 95% CL. The notation q 1 illustrates that only processes involving first generation 
quarks are considered. 
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